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We study the properties of rho mesons in nuclear matter by means of QCD sum 
rules at finite density. For increased sensitivity, we subtract out the vacuum contri- 
butions. With the spectral function as estimated in the literature, these subtracted 
sum rules are found to be not well satisfied. We suppose that Landau singularities 
from higher resonance states in the nearby region in this channel are the cause for 
this failure. 

1 Introduction 

It is believed that the mass and width of particles in a medium differ from their 
values in the vacuum. Such an effect might be observable in various physical 
systems such as the transition from a quark-gluon plasma to the hadronic phase 
in the early universe or in the interior of neutron stars. Furthermore, such 
conditions are created in heavy ion collisions at present and future colliders. 

In order to analyze these systems, it is necessary to obtain a good theoreti- 
cal understanding of the properties of the particles in a medium, for instance of 
the p-meson. There are various approaches and models in the literature which 
predict the behavior of the mass and width of the p-meson as function of the 
temperature and the density. We will use here the method o£QCD sum rules 
which has been quite successfully applied in the vacuum casecl. This approach 
was later extended to finite temperature and densitya. 

Since the medium breaks Lorentz invariance, more operators and unknown 
condensates enter on the QCD side of the sum rules (operator product expan- 
sion) and also the hadronic side (spectral function) is only poorly known in 
the medium. Nevertheless, one can perhaps gain insight into the chiral phase 
transition from experimental data on p-mesons in the medium, since these sum 
rules relate the mass and width of the p- meson to condensates like {qq} , which 
serves as an order parameter of chiral symmetry breaking. 

•Presented at the Conference on Strong and Electroweak Matter (SEWM 2000), Marseille, 
France, 14-17 June 2000. To be published in the Proceedings. Based on work done in 
collaboration with S. Mallik. 
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Most previous studies 13 found at zero temperature a drop of the p- meson 
mass of about 10-20% at nuclear saturation density n s , although in Rcf. □ no 
significant shift was obserised. In view of these conflicting results we want to 
scrutinize here and in Ref.u on some aspects of QCD sum rules in the medium. 

2 Operator product expansion and spectral representation 

We consider the ensemble average of the T-product of two vector currents 

7#(g) = i J d'xe^ (T (t£ (x)V v h '(0))) , (1) 

with V®(x) = q(x)j^(T a /2)q(x). The ensemble average is defined by (O) — 
Tv[exp(-(3{H - fi B N))0]/Tr[exp(-P(H - fJ, B N))] where (3 = 1/T. The nu- 
cleon chemical potential has been denoted by fi B and the QCD Hamiltonian 
and the nucleon number operator by H and N, respectively. The invariant 
decomposition of Tff is given by 

r^{q)=6 ab {Q liV T l + P^ v T t ), (2) 

with P M „ = -g^ v + {q^/q 2 ) - (q 2 / 'tf )u,j,u v and Q^ v = {q i / q 2 )^^, where 
Ufi = — ujq^/q 2 and q = \J uj 2 — q 2 . The amplitudes T/ it are functions of 
q 2 and u> — u ■ q. We have introduced the four- velocity of the medium u M , in 
order to formally restore Lorentz invariance which is broken by the medium. 

The sum rules are obtained by equating the short-distance expansion of the 
product of the two currents in Eq. ([!]), which results in a series of condensates of 
QCD operators, with the spectral representation for the invariant amplitudes 
Tij, where various intermediate hadronic states will contribute. 

The short distance expansion of 7T™ including all operators up to dimen- 
sion four can be found in Ref.u. The availability of the four- vector allows 
one to construct new operators as compared to the vacuum case. In addition 
to the usual operators 1, qq — uu + dd, G 2 = ^-G"G^ va , there are two 
new operators G« = u^O^u" and @ g = u^Q^t? . Here G° v , a = 1,...,8 
are the gluon field strengths and a s — g 2 /4tt, g being the QCD coupling con- 
stant. 0*$ are the quark and the gluon parts of the traceless stress tensor 
®U = <jilf*D u q - fg^qq and 0£„ = -G^G u Xc + \g^G c afi G a ^, where m is 
the quark mass in the limit of SU(2) symmetry. We will take into account the 
mixing of the operators Q q and 9 under the renormalization group. 

At finite temperature and chemical potential, we can use the Landau rep- 
resentational of the amplitudes T/ jt , which is a spectral decomposition in qo at 
fixed \q\ 

TiAlo, m = - l +C ° ^ ImrLt(g °' lg1) tanh(/3^/2), (3) 
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up to subtraction terms. There will be contributions in the integral from 
various intermediate hadronic states, first of all from the p-meson. The effects 
of the medium can be parametrized by employing the operator relation 

VZ{x)=m*F* p p%{x), (4) 

where m* and F* denote the in-medium mass and width of the p-meson. This 
will generate the usual (5-function contribution in ImTi t- 

Moreover, in the nuclear medium there will be a contribution from NN 
intermediate states. In the vacuum this contribution is small, coming from the 
cut beginning at threshold, q$ = Am? N + \q\ 2 . However, in the nuclear medium 
the currents can also interact with real nucleons to give rise to a short cut 
around the origin, — \q\ < qa < The evaluation of this contribution can 

be found in Ref.LJ. 

Following Refs.i'§ we ignore the states NN* with resonances N* . Fur- 
thermore, below we will restrict ourselves to sum rules at finite nuclear density 
but at zero temperature, therefore the 2tt contribution will be eliminated after 
subtracting out the corresponding vacuum sum rules. 

3 Sum rules 

We obtain the sum rules by equating the spectral representation and the opera- 
tor product expansion for the two invariant amplitudes T/. f . As usually done in 
the literature^, we take the Borel transform of both sides in order to enhance 
the contribution from the lowest lying resonance, here the p-meson, and to 
suppress the contributions from higher dimensional operators. In general, this 
can only be achieved within a xiertain region for the so called Borel parameter 
M . In contrast to earlier workHu we subtract the vacuum sum rules, assuming 
that the contribution from the QCD continuum will practically drop out in 
this way. For T — > and fig > the 2tt contribution will also cancel and 
we obtain the following sum rules in the limit \q\ — > where the expressions 
simplify considerably: 

I ^ 1 f iA , , --l, / 4mi, 2m%. JO) , . 
Ffe * + — 2 Ml -e^l- + = L| , (5) 

1 fiA , / A m 2 2TT7 2 

mfFfe ^_ 1 / dsse -^l 1 _^N {1 + ^N ) = _ {0)A6) 
J 4m 2 N V S S 

with 

(O) = \m{qq) + ^ + A ({Q) + A(M 2 )(^6^ - &>?) , (7) 
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where (O) = (O) — (0|O|0) and 9 = Q q + Q 9 . The mixing under the renormal- 
ization group is taken into account by \(M 2 ) = (a s ([i 2 ) /a s (M 2 ))~ d / 2b with 
d = f + ?!/) and 6 = 11 - §n/. We will take n = 1 GeV below. 

In the linear density approximation, we expand the mass and width and 
the condensates up to first order in the nucleon number density n. 
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m p (l + o— ), F* = F p (l + &— ), {0) = Cn, 
n s 1 n s 



(8) 



where n s = (110 MeV) 3 denotes the nuclear saturation density. The coefficient 
C in Eq. (0) is given by 



n a 1 ( 
C=---{m N 



a) + -m N 



A q + A 9 + X(M 2 )(-A q - A 9 ) 



(9) 



where a — (p\m(uu + dd)\p) /(2mjv) — 45 MeV denotes the sigma termEl. The 
constants A 9,9 are defined by {p\Q q ^ 9 \p) — 2A q > 9 {p^p v — jg^p 2 )- From parton 
distribution functions □ one obtains the values A q = 0.62 and A 9 = 0.35. 
From the sum rules (|5|) and (||) we finally obtain the expressions 
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4 Discussion and Conclusions 

In Fig. [I] we have plotted the coefficients a and b as function of the Borel 
parameter M. We note that there is no sign of constancy of a and b in any 
region of M. We thus conclude that the sum rules in Eqs. (|5|) and (||) cannot 
give any reliable information about the density dependence of m p and F p . 

We have included only operators up to dimension four in the sum rules 
above, but taken into account their mixing which is numerically relevant for 
M 2 ^ /j, 2 . Contributions from higher dimensional operators should, however, 
be relatively small for M 2 > 1 GeV 2 . A more detailed analysis is in progress^. 

Therefore, the failure of our subtracted sum rules can presumably be traced 
to the hadronic spectral side which is not adequately saturated. In fact, higher 
resonance states NN* will also contribute a cut for q 2 < (m* N — win) 2 . The 
problem with the inclusion of such contributions is, however, that more un- 
known couplings and masses will enter in the medium. 

The sum rules at finite density are given by the vacuum sum rules per- 
turbed by small terms proportional to the density. Since the vacuum sum rules 
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Figure 1: Coefficients of n/n s in the expansion of m* and F* from Eq. ^ as function of 
the Borel parameter M. 

are stableS, this guarantees-.the stability of the sum rules at finite density as 
observed in previous workcKa. In contrast, the subtracted sum rules above are 
much more sensitive to errors in or omissions of any terms. 
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